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Introduction 

The  United  States  Army  and  Air  Force  use  MIL-T-83133C  grade  JP8  as  the  main  source  of 
energy  for  military  land  vehicles  and  aircraft.  JP-8  is  very  similar  to  Jet  A-l  with  three 
military  specified  additives,  fuel  system  inhibitor,  corrosion  inhibitor/lubricity  improver 
and  static  dissipater.  The  conversion  to  JP-8  occurred  primarily  to  improve  the  safety  of 
aircraft,  although  the  “single  fuel  for  the  battlefield”  concept  is  centered  on  the  use  of 
aviation  kerosene  in  all  Air  Force  and  Army  aircraft  and  ground  vehicle. 

Detailed  chemical  kinetic  mechanisms  that  describe  combustion  of  many  of  the 
components  in  JP-8  are  not  available  and  are  unlikely  to  be  developed  in  the  near  future. 
Hence  there  is  a  need  to  study  the  characteristics  of  JP-8  experimentally.  One  of  the 
important  thermo-physical  properties  of  JP-8/air  mixtures  is  its  laminar  flame  speed.  There 
are  many  more  difficulties  in  the  accurate  measurement  of  the  laminar  burning  speeds  of 
liquid  fuels  such  as  JP-8  than  that  of  gaseous  fuels.  The  process  of  creating  a  precise  initial 
mixture  of  the  fuel/oxidizer/diluents  can  present  many  challenges.  In  this  research  the 
unique  experimental  facility  in  Northeastern  University  allows  to  prepare  correct  mixture  at 
high  temperatures  and  pressures  to  perform  the  experiments. 

The  following  is  the  final  report  on  our  developments  and  accomplishments  for  the 
three  years  period.  In  this  report,  a  brief  description  of  the  experimental  facilities  which 
includes  two  combustion  chambers,  spherical  and  cylindrical,  optical  set  up,  a  high 


temperature  oven  will  be  presented.  The  structure  of  JP8-C>2-diluent  flames  have  been 
studied,  laminar  burning  speeds  of  JP8-C>2-He,  JP8-C>2-Ar  and  JP8-air  mixtures  have  been 
calculated  and  a  power  law  fit  has  been  developed  for  the  calculation  of  laminar  burning 
speed  of  these  fuels  at  different  conditions.  Also,  the  cell  formation  and  instability  of  these 
mixtures  have  been  investigated. 

Autoignition  and  explosion  limits  of  JP-8/air  mixtures  have  been  measured  by 
performing  some  experiments  in  the  spherical  vessel.  The  experiments  were  perfonned  at 
high  initial  pressures  and  temperatures  to  observe  autoignition.  The  explosion  limits 
temperatures  and  pressures  were  measured  at  different  equivalence  ratios. 

Experimental  Set  up 

Burning  speed  measurements  were  made  in  the  existing  spherical  combustion 
chamber.  The  spherical  chamber  consists  of  two  hemispheric  heads  bolted  together  to  make 
a  15.4  cm  inner  diameter  sphere.  The  chamber  was  designed  to  withstand  pressures  up  to 
425  atm  and  is  fitted  with  ports  for  spark  electrodes,  diagnostic  probes,  and  ports  for  filling 
and  evacuating  it.  A  thennocouple  inserted  through  in  one  of  the  chamber  ports  was  used  to 
check  the  initial  temperature  of  the  gas  inside  the  chamber.  A  Kistler  603B 1  piezo-electric 
pressure  transducer  with  a  Kistler  501  OB  charge  amplifier  was  used  to  obtain  dynamic 
pressure  vs.  time  records  from  which  the  burning  speed  was  detennined.  Ionization  probes 
mounted  flush  with  the  wall  located  at  the  top  and  bottom  of  the  chambers  were  used  to 
measure  the  arrival  time  of  the  flame  at  the  wall  and  to  check  for  spherical  symmetry  and 
buoyant  rise. 

Spherical  vessel  is  housed  in  an  oven  which  can  be  heated  up  to  500  K.  Liquid  fuel 
is  stored  in  a  1 15  cc  heated  vessel  and  is  transferred  through  a  heated  line  inside  the  oven  to 
the  spherical  chamber.  Several  thermocouples  are  located  on  the  line  from  the  fuel 
reservoir  to  the  vessel  to  monitor  temperature  of  the  fuel  passageway.  A  heated  strain 
gauge  (Kulite  XTE-190)  in  the  oven  is  used  to  measure  partial  pressure  of  fuel  in  the 
vessel. 

The  companion  cylindrical  chamber  is  made  of  SAE4140  steel  with  an  inner 
diameter  and  length  of  133.35  mm.  The  two  end  windows  are  34.93  mm  thick  Pyrex  with  a 
high  durability  against  pressure  and  temperature  shocks  as  well  as  having  very  good  optical 


4 


properties.  A  Z-type  Schlieren/Shadowgraph  ensemble  has  been  set  up  to  visualize  the 
flame  propagation.  A  high  speed  CMOS  camera  (MotionXtra  HG-LE,  Redlake  Inc.)  with  a 
capture  rate  of  up  to  40,000  frames  per  second  is  placed  very  close  to  the  focal  point  of  the 
second  mirror.  Two  band  heaters  and  a  rope  heater  wrapped  around  the  cylindrical  vessel 
are  used  to  heat  up  the  inside  temperature  of  the  vessel  up  to  500  K.  This  chamber  is 
equipped  with  a  heated  liquid  fuel  line  system,  a  pressure  strain  gauge  and  thermocouples 
similar  to  the  spherical  vessel.  The  oven  was  omitted  to  permit  flame  observation  for  this 
application. 

Gas  handling  system  used  with  these  facilities  consists  of  a  vacuum  pump  for 
evacuating  the  system  and  a  valve  manifold  connected  to  gas  cylinders  for  preparation  of 
the  fuel/oxidizer/diluent  mixtures.  Partial  pressures  of  the  fuel  mixtures  were  measured 
using  Kulite  strain  gauge  pressure  transducers  in  the  0-15  psi  range.  Two  spark  plugs  with 
extended  electrodes  were  used  to  ignite  the  mixture  at  the  desired  location  in  the  chambers. 
An  electronic  ignition  system  controlled  by  the  data  acquisition  program  provides  a  spark 
with  the  necessary  energy.  The  data  acquisition  program  synchronizes  the  ignition  with  the 
dynamic  pressure  recording  and  Schlieren/Shadowgraph  photography. 

The  data  acquisition  system  consists  of  a  Data  Translations  16  bit  data  acquisition 
card,  which  records  the  pressure  change  of  the  combustion  event  at  a  rate  of  250  kHz.  The 
analog  to  digital  converter  card  receives  the  pressure  signal  from  the  charge  amplifier  and 
the  signals  from  the  ionization  probes.  All  signals  are  recorded  by  a  personal  computer  and 
an  output  data  file  is  automatically  generated.  The  output  data  files  include  the  dynamic 
pressure  and  its  corresponding  time. 

Test  procedure  begins  by  evacuating  the  vessel  and  gas  handling  system  using  the 
vacuum  pump.  The  chamber  then  is  filled  with  JP8  vapor  to  the  desired  pressure  and  air  is 
added.  The  vessel  and  the  fuel  tank  are  at  the  same  temperature  during  the  filling.  After  the 
chamber  is  filled  with  the  proper  mixture,  several  minutes  are  allowed  for  the  system  to 
become  quiescent  before  it  is  ignited.  This  will  prevent  any  turbulence  inside  the  vessel. 
Six  thermocouples  on  the  liquid  line  are  used  to  make  sure  that  temperature  along  the 
filling  line  is  never  below  the  condensation  temperature  for  JP8.  At  least  three  runs  at  each 
initial  condition  were  made  to  provide  a  good  statistical  sample.  Based  on  statistical 
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analysis,  it  was  found  that  three  runs  are  sufficient  to  achieve  a  95%  confidence  level. 
Figure  1  shows  the  schematic  of  two  combustion  chambers. 


Theoretical  Model 

The  theoretical  model  used  to  calculate  the  burning  speed  from  the  pressure  rise  is  based  on 
one  previously  developed  by  Metghalchi  and  co-workers  [1-5]  and  has  been  modified  to 
include  corrections  for  heat  losses  to  electrodes  and  radiation  from  the  burned  gas  to  the 
wall  as  well  as  the  temperature  gradient  in  the  preheat  zone.  It  is  assumed  that  gases  in  the 
combustion  chamber  can  be  divided  into  burned  and  unbumed  regions  separated  by  a 
reaction  layer  of  zero  thickness  as  shown  schematically  in  Figure  2.  Burned  gas  in  the 
center  of  chamber  is  divided  into  n  number  of  shells  where  the  number  of  shells  is 
proportional  to  the  combustion  duration.  Burned  gas  temperatures  in  shells  are  different 
from  each  other  but  burned  gases  are  in  chemical  equilibrium  in  each  shell  at  the  given 
temperature.  The  burned  gas  is  surrounded  by  the  unburned  gas  in  preheat  zone  ( 5  ,  )  with 

variable  temperature.  Core  unburned  gas  with  uniform  temperature  surrounds  the  preheat 
zone  gases.  A  boundary  layer  ( 8hl )  separates  core  unburned  gas  from  vessel  wall.  It  is 

further  assumed  that:  the  burned  and  unbumed  gases  are  ideal,  the  unbumed  gas 
composition  is  frozen,  the  pressure  throughout  the  chamber  is  uniform,  compression  of 
both  burned  and  unbumed  gases  is  isentropic,  and  the  heat  flux  due  to  the  temperature 
gradient  in  the  burned  gas  is  negligible. 

Radiation  in  the  participating  media  is  a  well  known  area  and  has  received  a  lot  of  attention 
in  the  past  fifty  years.  The  ideal  approach  to  calculate  radiative  energy  transfer  through  the 
hot  gases  in  an  enclosure  is  to  use  the  relations  for  absorption,  emission  and  scattering 
coefficients  and  their  corresponding  configuration  factors.  However,  simple  solutions  for 
temperatures  and  heat  flows  are  usually  not  possible  for  equations  of  radiative  transfer 
combined  with  energy  conservation  relations  including  the  correlations  of  participating 
media.  For  some  conditions  however,  simplification  can  be  made  that  will  provide  useful 
results.  Among  all  different  approaches,  optically  thin  and  thick  approximation  which  their 
limits  are  examined  in  [6-7]  has  been  predominantly  used. 

In  our  study  it  has  been  assumed  that  gases  in  the  chamber  are  optically  thin.  This 
assumption  is  valid  for  the  flames  under  study  since  even  the  species  with  largest 
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absorption  coefficients  do  not  play  a  major  role  in  reabsorption.  The  reason  is  that,  the 
product  of  the  largest  attenuation  coefficient  by  the  dimension  of  the  combustion  bomb  is 
much  less  than  unity.  This  argument  has  been  also  supported  by  Vranos  and  Hall  [8]  who 
has  done  a  theoretical  analysis  on  the  effect  of  thennal  radiation.  They  used  gas  band 
radiation  description  and  allowed  for  reabsorption.  Based  on  their  detailed  calculation, 
reabsorption  becomes  important  only  at  very  low  strain  rates  and  the  correction  for  strain 
rates  of  30  s'1  and  below  is  less  than  0.2%.  The  effect  of  radiation  on  temperature  becomes 
more  pronounced  as  strain  rate  decreases.  In  our  experiments  with  different  fuel/oxidizer 
mixtures  strain  rate  at  its  lowest  value  is  not  less  than  30  s'1  and  hence  the  correction  for 
reabsorption  is  very  insignificant. 

It  is  assumed  that  there  is  an  energy  transfer  process  from  the  mixture  to  the  spark 
electrodes.  Different  spark  electrodes  have  been  tested  with  various  thicknesses  to 
investigate  the  effect  of  energy  transfer  on  the  burning  speed  but  it  was  observed  that  it  did 
not  affect  significantly.  However,  the  energy  transfer  term  to  the  electrodes  is  involved  in 
the  model.  The  thinnest  electrodes  have  been  used  to  reduce  energy  transfer. 

For  the  conditions  of  interest  in  the  present  work,  all  these  assumptions  have  been  validated 
by  numerous  experiments  in  constant  volume  chambers  and  internal  combustion  engines 
carried  out  over  the  past  several  decades. 

For  spherical  flames,  the  temperature  distribution  of  the  gases  in  the  combustion  chamber 
and  the  burned  gas  mass  fraction  can  be  determined  from  the  measured  pressure  using  the 
equations  for  conservation  of  volume  and  energy  together  with  the  ideal  gas  equation  of 
state 

pv  =  RT  (1) 

where  p  is  the  pressure,  v  is  the  specific  volume,  R  is  the  specific  gas  constant  and  T  is  the 
temperature. 

The  mass  conservation  equation  is 

m  =  mh  +  mu  =  pbVb  +  puVu  =  pt  (Vc-Ve)/RTt  (2) 

where  m  is  the  mass  of  gas  in  the  combustion  chamber,  mb  is  the  burned  gas  mass,  mu  is  the 
unbumed  gas  mass,  Vc  is  the  volume  of  the  combustion  chamber,  Ve  is  the  electrode 
volume,  and  the  subscript  i  denotes  initial  conditions.  Subscripts  b  and  u  represent  burned 
and  unburned  conditions  respectively,  p  is  average  density  and  V  is  the  volume  of  the  gas. 
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The  total  volume  of  the  gas  in  the  combustion  chamber  is 

Vi=Vc-Ve=Vh+Vll 

The  energy  conservation  equation  is 

E.-Q  -Q  -Q  =  E,  +  E 

i  z^e  w  b  u 


(3) 

(4) 


where  E,  is  the  initial  energy  of  the  unbumed  gas,  Qw  is  the  conduction  heat  loss  to  the 
wall,  Qe  is  the  conduction  heat  loss  to  the  electrodes,  Q,  is  the  heat  loss  due  to  radiation 
from  the  burned  gas. 

In  the  case  of  rapidly  increasing  pressure  such  as  that  occurring  during  constant  volume 
combustion,  the  terms  representing  compression  work  on  the  boundary  layer  may  be 
neglected  and  the  resulting  equations  are 

Q,=pV,bKn-i)  =  E,l  (5) 

Q,,=pVwtHr,-i)  =  E,„  (0 


The  radiation  heat  loss  from  the  burned  gas  was  calculated  using 

Q,=iQM')dt'  =  4arViaTl4  (7) 

where  ap  is  the  Planck  mean  absorption  coefficient  and  a  is  the  Stefan-  Boltzman  constant. 
Finally,  we  solve  the  following  equations  of  volume  and  energy  simultaneously 

£6(vfa<T>)-vJ^'  =  v,. -vus  +(Veh  +Vwb  +  Vph)/m 

=e<  ~eus  +(pVp'r/^u  -b-QrV™  (9) 

where  Vi  =  (Vc-Ve)/m  and  c,  =  E,/m  are  the  initial  specific  volume  and  energy  of  the 
unburned  gas  in  the  chamber. 

Equations  8  and  9  contain  the  three  unknowns:  p,  Xb(p),  and  Tb(r,t).  Using  measured,  p(t), 
as  a  function  of  time,  they  can  be  solved  numerically  using  the  method  of  shells  to  obtain 
the  burned  mass  fraction,  Xb(t),  as  a  function  of  time  and  the  radial  temperature  distribution 
T(r,t).  The  mass  burning  rate,  mh  =  mxh  can  be  obtained  by  numerical  differentiation  of 


Xb(t).  The  thermodynamic  properties  of  the  burned  and  unburned  used  in  the  calculations 
were  obtained  from  the  JANAF  Tables  and  STANJAN  code. 


If  the  specific  heats  of  the  burned  and  unburned  gases  are  assumed  constant  we  can 
simplify  Eq.  9 

Pivjin-D-  v.  Hr,  -D)-  K„  +  + Ar  -Q»m  (10) 

(Yu  ~  7)  ( Yu  - 

where  hub  =  hfu  -  ha  is  the  difference  between  the  enthalpies  of  formation  of  the  unburned 
and  burned  gas  at  zero  degrees  Kelvin.  Equations  8  and  10  are  identical  to  those  previously 
reported  in  last  year’s  reports  with  correction  for  radiation  heat  loss,  conduction  heat  loss  to 
the  electrodes,  and  the  displacement  thickness  of  the  preheat  zone  added.  Given  the 
pressure,  they  can  be  solved  analytically  to  obtain  an  excellent  approximation  for  Xb  and 
(RT)bs  =  pvbs.  However,  they  do  not  give  the  temperature  distribution  in  the  burned  gas.  It 
is  of  interest  to  note  that  if  the  small  correction  terms  for  the  preheat  zone  and  the  heat  loss 
to  the  electrode  are  neglected,  xb  and  (RT)bs  =  pvbs  depend  only  on  the  pressure  and  are 
independent  of  the  flame  shape. 

Burning  speed  may  then  be  defined 

Sh  =  mb/  puAb  =  mxl  puAb  (11) 

where  A/,  is  the  area  of  a  sphere  having  a  volume  equal  to  that  of  the  burned  gas.  This 
expression  is  valid  for  smooth,  cracked,  or  wrinkled  flames  of  any  shape. 

Experimental  Results 

A  set  of  experiments  were  performed  to  study  the  structure  of  JP8-02-diluent  premixed 
flames  and  ultimately  calculate  the  burning  speed  of  the  laminar  cases.  Power  law  fits  were 
developed  for  the  laminar  burning  speeds  by  least  square  method.  The  performed 
experiments  are: 

1-  JP8-Air  experiments  in  equivalence  ratio  range  of  0.8- 1.2  and  initial  pressure  range 
of  1-5  atm 

1-1:  Burning  speed  is  calculated  and  power  law  curve  fit  is  developed  by  least 
square  method  for  laminar  cases.  The  laminar  burning  speed  fit  are  for  the  range  of 
0.8<  q><l  .2,  500<T<675  K  and  1<P<7  atm. 

1-2:  Pictures  have  been  taken  in  those  conditions. 

2-  JP8-OyHe  experiments  in  equivalence  ratio  range  of  0.8-1  and  initial  pressure 
range  of  1-7  atm 
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2-1:  Burning  speed  is  calculated  both  in  laminar  and  non-laminar  cases.  Power  law 
fit  is  developed  for  laminar  for  0.8<cp<  1 ,  500<T<800  K,  and  1<P<20  atm. 

2- 2:  by  shadowgraph  system  a  set  of  pictures  have  been  taken  in  pressure  range  of 
1-30  atm  and  equivalence  ratio  range  of  0.8-1. 

3-  JP8 -Qi-Ar  experiments'. 

3- 1:  set  of  flame  propagation  pictures  are  in  range  of  0.8- 1.2  and  pressure  range  of 
1-23  atm. 

4^  Instability  analysis : 

Thermo-diffusive  and  hydrodynamic  instability  and  cell  fonnation  analysis  have 
been  done  for  JP8  at  different  initial  pressures,  diluents  and  equivalence  ratios. 
Some  important  properties  of  the  flame  have  been  calculated  (flame  thickness, 
thennal  diffusion  coefficient).  The  analyses  show  systematic  and  reproducible 
results. 

5z  Diesel  Fuel  experiment: 

Burning  speed  of  diesel-air  mixtures  is  calculated  for:  0.8<cp<l ,  500<T<760K  and 
1<P<6  atm.  Flame  pictures  for  atmospheric  initial  pressure  have  been  taken. 

6z_  JP-8/air  autoignition  experiments'. 

The  autoignition  and  explosion  limits  of  JP-8/air  mixtures  were  tested  by 
conducting  high  temperature  and  high  pressure  experiments  in  spherical  vessel.  The 
autoignition  temperature  and  pressure  of  theses  mixture  were  measured  in  various 
equivalence  ratios. 


JP8-Air  experiments 

JP8-air  experiments  have  been  performed  to  measure  the  burning  speed  and  to  study  flame 
structures  of  those  mixtures.  Figure  3  shows  the  flame  propagation  of  JP8-air  mixtures  at 
cp=  0.8  and  three  different  initial  pressures  which  are  P;=l,  2.5  and  5  atm.  It  is  clear  that  as 
the  pressure  of  the  chamber  increases  the  cells  are  formed.  Figure  4  presents  the  flame 
propagation  of  JP8-air  at  cp=  0.9  at  the  same  initial  pressures.  Compared  with  <p=  0.8 
flames,  the  cells  are  formed  earlier.  The  same  result  is  obtained  by  looking  the  figures  5 
and  6  for  equivalence  ratios  of  1  and  1.2.  Two  major  results  are  concluded  in  the  pictures: 
instability,  and  ultimately  turbulence,  originates  from  cells  in  the  flame  which  fonn  with 
increasing  pressure;  cell  formation  develops  earlier  as  the  mixture  gets  richer. 
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Table  1  presents  the  cell  formation  pressures  and  temperatures.  The  critical  pressure  has 
been  measured  by  the  pressure  transducer,  and  the  critical  temperature  has  been  calculated 
by  isentropic  compression  assumption  of  unburned  gas  zone. 


(12) 


where  Tcr  and  Pcr  are  cell  fonnation  temperature  and  pressure  respectively  and  Pi  and  T,  are 
initial  pressure  and  temperature  respectively,  y  is  the  ratio  of  constant  pressure  heat 
capacity  to  the  constant  volume  heat  capacity. 

The  same  experiments  were  used  to  calculate  the  burning  speed  of  JP8-air  mixtures.  Figure 
7  show  a  typical  pressure  curve  versus  time.  The  black  solid  curve  is  the  pressure  and  the 
colored  lines  are  for  the  ionization  probes  sample.  Ionization  probes  are  located  on 
spherical  vessel  walls  to  pinpoint  the  arrival  of  the  flame.  The  discontinuous  part  of 
colored  lines  show  the  moment  when  the  flame  has  touched  the  wall,  if  the  discontinuity 
be  observed  at  the  same  time,  it  is  understood  that  the  flame  is  symmetrical  and  buoyancy 
effects  are  negligible. 

The  laminar  burning  speed  has  been  calculated  and  Figure  8  presents  the  laminar  burning 
speed  curves  of  the  JP8-air  mixtures  at  initial  pressure  of  1  atm.  The  burning  speed  data 
are  plotted  on  isentropic  lines  versus  unbumed  gas  temperature.  As  it  is  expected  in  lean 
conditions  the  burning  speed  increases  by  increasing  the  equivalence  ratio.  Figure  9  and  10 
show  the  laminar  burning  speed  curves  of  the  JP8-air  mixtures  at  initial  pressure  of  2  and  5 
atms.  The  same  trend  is  observed  but  by  increasing  the  pressure  the  burning  speed  of  the 
mixture  decreases.  A  power  law  fit  has  been  developed  by  least  square  method  to  calculate 
the  burning  speed  of  JP8-air  mixtures.  The  form  of  equation  is: 

~\a  r  ~\P 

Sb  =  Sb0(l  +  a](l-<f>)  +  a2(l-<f>)2)  ^  (13) 

_  u0  J  L  -*  0  _ 


where  Sbo  is  the  burning  speed  at  reference  point  (P0  =  1  atm,  Tuo  =  450  K  and  <j>  =  1)  in 
cm/s,  (f>  is  the  mixture  equivalence  ratio,  Tu  is  the  unburned  gas  temperature  in  K,  Tu0  is 
the  reference  temperature  and  equal  to  450  K,  P  is  the  mixture  pressure  in  atm  and  Po  is  the 
reference  pressure  and  equal  to  1  atm.  The  fitted  constants,  ai,  a2,  a  and  //  are  shown  in 
table  2. 
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JP8-02-Ar  experiments 

By  looking  at  figures  3-6  it  is  understood  that  cells  are  formed  in  JP8-air  mixtures 
early  and  due  to  this  fact,  developed  power  law  fit  for  laminar  flame  can  not  cover  a  wide 
range  of  temperatures  and  pressures.  It  was  planned  to  replace  the  diluent  of  the  original 
mixture  (Nitrogen)  with  Argon  to  see  the  effect.  Figure  1 1  shows  the  comparison  between 
the  flame  images  of  JP8-air  and  JP8-02-Ar  mixtures.  As  it  is  seen  in  the  pictures,  Argon  is 
not  a  stabilizing  factor,  causing  cell  formation  earlier  than  Nitrogen 

JP8-02-He  experiments 

The  diluent  was  replaced  with  helium  to  see  its  effect  on  the  structure  of  the  flame.  Figure 
12  shows  the  JP8-02-He  flames  at  three  different  equivalence  ratios  0.8,  0.9,  and  1  initial 
temperature  of  463  K  and  initial  pressure  of  4  atm.  The  pictures  demonstrate  that  helium  as 
diluent  can  be  a  stabilizing  factor  compared  with  Nitrogen  and  Argon.  It  can  be  due  to 
larger  Lewis  number  and  smaller  expansion  ratio  of  JPS-CL-Hc  mixtures.  Table  3  presents 
the  pressure  and  temperatures  in  which  the  cells  were  fonned.  The  data  show  that  the  cells 
are  formed  at  a  specific  pressure  which  is  a  function  of  equivalence  ratio  and  diluent  type. 
The  cell  fonnation  temperature  can  be  calculated  by  the  equation  (12).  The  same 
experiments  were  done  in  the  spherical  vessel  to  calculate  the  laminar  burning  speed  of 
JP8-02-He.  Figure  13  shows  the  laminar  burning  speeds  of  JP8-02-He  mixtures  at  initial 
pressure  of  1  atm  versus  unburned  gas  temperatures.  As  it  is  seen  in  the  lean  mixtures, 
laminar  burning  speed  increases  as  the  mixture  becomes  richer.  Figure  14  presents  the 
laminar  burning  speed  at  initial  pressure  of  2  atm,  the  same  trend  is  observed  without  any 
significant  difference  from  Figure  13.  Figures  15,  16  and  17  are  the  laminar  burning  speed 
curves  at  initial  pressures  of  4,  5,  7  atms,  respectively.  Contrary  to  JP8-air  mixtures,  initial 
pressure  does  not  play  any  role  on  the  laminar  burning  speed  of  the  mixture.  A  polynomial 
function  has  been  fitted  by  least  square  method  for  the  laminar  burning  speed  of  JP8-02-He 
mixtures.  The  form  of  the  function  is  the  same  as  equation  (13).  The  corresponding 
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coefficients  can  be  found  at  Table  4.  As  it  can  be  seen  the  pressure  exponent  is  almost  zero 
making  burning  speed  independent  of  pressure. 


Instability  analysis 

A  set  of  experiments  were  performed  on  JP8-C>2-diluent  mixtures  to  detect  the  onset 
of  instability  and  effective  parameters  in  more  detail.  The  strategy  is  to  investigate  the 
occurrence  of  instability  by  performing  many  experiments.  The  changing  parameters  are 
equivalence  ratio,  initial  pressure  and  diluent  type.  The  criterion  for  the  instability  is  the 
formation  and  evaluation  of  small  cells  which  develop  by  the  propagation  of  the  flame. 
Critical  pressures  and  temperatures  that  cells  are  formed  are  tabulated  in  tables  5  and  6.  It 
can  be  seen  that  the  important  factor  in  the  cell  formation  and  instability  of  the  JP8-O2- 
diluent  flames  are  pressure  and  flame  thickness.  An  estimation  of  flame  thickness  can  be 
calculated  by  the  following  formula: 

5«—  (14) 

Su 

where  8 ,  a  and  Su  are  flame  thickness,  thermal  diffusivity  and  the  burning  speed  of  the 
mixture  respectively.  As  it  is  seen,  the  cell  formation  pressure  is  independent  of  initial 
pressure  and  the  instability  temperature  just  depends  on  the  pressure  due  to  isentropic 
compression  assumption.  Another  result  is  that  in  lean  mixtures  flame  becomes  more 
stable  and  the  cell  fonnation  pressure  increases.  There  is  not  any  huge  change  in  the  flame 
thickness,  but  by  non-dimensionalization  of  the  flame  thickness  versus  cell  formation 
radius  (Rcr),  a  systematic  trend  is  observed  which  is  the  increasing  S/Rcr  by  the  increasing 

of  the  equivalence  ratio  (8/Rcr  is  the  Peclet  number). 

The  effect  of  diluent  was  investigated  next.  The  diluent  (N2)  was  replaced  with 
Helium  to  see  the  effect.  It  was  observed  that  in  the  mixtures  with  helium  as  diluent,  the 
cell  fonnation  was  delayed  because  of  higher  instability  pressure.  Table  6  presents  the 
results  of  these  experiments. 

The  difference  between  mixtures  with  different  diluents  can  be  attributed  to  the  heat 
conductivity.  One  of  the  most  important  factors  in  flame  stabilization  is  the  energy  and 
mass  transfer  properties  of  the  flame.  Energy  transfer  is  characterized  by  heat  conduction 
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and  mass  transfer  is  specified  by  Dm,  the  mass  diffusion  coefficient.  Combining  both  of 
these  parameters  produces  the  Lewis  number: 


Le  = 


k 

PCPD, 


(15) 


where  k  is  the  heat  conductivity  of  the  mixture,  p  is  the  density,  cp  is  the  heat  capacity,  and 
Dm  is  the  mass  diffusion  coefficient  of  deficient  reactant  in  the  mixture. 

It  is  well  established  that  mixtures  with  higher  Lewis  numbers  are  more  thermo- 
diffusibly  stable  [9-12].  A  FORTRAN  program  is  developed  which  is  joined  with 
CHEMKIN  to  calculate  the  Lewis  number  of  different  mixtures.  As  there  is  not  any 
released  transport  properties  data  for  JP8  (which  is  needed  as  input  for  CHEMKIN),  Lewis 
numbers  of  propane  (C3H8)  were  calculated  with  different  diluents  at  equal  temperatures 
and  pressures.  It  was  seen  that  Lec  H%_He  >  Lec  Hg  Ar  >  Lec  H  Ni  which  roughly  confirms 

the  stability  of  JP8-02-diluent  mixtures.  In  addition,  flame  thickness  of  CfitL/CF/diluent 
flames  has  been  calculated  and  it  has  been  seen  thatbc  Hs_He  >  5C3Hs_Ni  >  8ch  Earlier 

cell  formation  of  JP8-Ar  flames  can  be  due  to  thinner  flame  thickness.  Figure  18  shows 
the  flame  propagation  pictures  of  JP8-O2  mixtures  with  three  different  diluents  and  the 
analyses  are  confirmed.  It  is  obvious  that  the  mixtures  with  Argon  as  diluent  are  the  most 
sensitive  mixtures  for  cell  fonnation  due  to  thinner  flame  thickness  compared  to  the  other 
ones.  The  mixtures  with  Helium  as  diluent  are  the  most  stable  ones  due  to  higher  Lewis 
number  and  thicker  flames. 

The  reason  for  the  mentioned  phenomena  is  complex  because  of  the  involvement  of 
both  hydrodynamics  and  thermo-kinetics.  The  next  task  of  this  research  is  the  calculation 
of  some  other  important  factors  like  expansion  ratio  and  Lewis  number  of  the  mixtures  to 
extract  a  rational  trend. 


Diesel  Fuel  Experiments 

Diesel-air  experiments  have  been  performed  at  atmospheric  initial  pressure.  Flame 
pictures  demonstrate  that  cells  are  formed  at  early  times  after  ignition.  It  shows  that  diesel- 
air  mixtures  are  more  sensitive  to  instabilities  and  cell  fonnation  than  other  kinds  of  flames. 
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The  reason  for  this  is  not  well  understood,  since  the  structure  of  diesel  fuel  is  not  defined 
yet.  Figure  19  shows  the  flame  propagation  of  diesel-air  mixtures  at  atmospheric  initial 
pressure,  stoichiometric  equivalence  ratio  and  initial  temperature  of  463  K.  It  is  obvious 
that  the  tolerance  of  diesel-air  mixtures  to  instability  is  less  than  JP8-air  mixtures.  The 
advantage  or  disadvantage  of  turbulence  occurrence  in  the  flames  depends  on  the 
application  which  is  expected  of  the  combustible  device.  Diesel  fuel  is  going  to  be 
replaced  with  JP8  in  some  of  machines.  With  respect  to  the  application  of  the  device,  this 
should  be  seen  whether  lower  cell  formation  tendency  of  JP8  is  a  positive  sign  or  not. 

The  same  experiments  were  performed  in  spherical  vessel  and  burning  speed  of  the  mixture 
was  calculated.  Figure  20  shows  the  burning  speeds  of  diesel/air  mixtures  at  three  different 
equivalence  ratios  and  initial  pressure  of  1  atm.  As  expected  this  is  seen  that  at  lean 
conditions,  by  increasing  the  equivalence  ratio  burning  speed  increases. 

JP-8/air  autoignition  experiments 

The  experiments  have  been  performed  at  four  initial  pressure  conditions  of  8,9,10,  and  11 
atm  respectively.  The  initial  temperature  was  constant  at  215  °C  to  ensure  that  it  meets  the 
boiling  temperature  of  JP-8  [13],  and  the  experiments  were  carried  out  using  three 
equivalence  ratios  of  0.85,1,  and  1.2. 

The  initial  conditions  are  very  important  since  these  experiments  are  strongly 
equivalence  ratio  and  temperature  dependent.  The  procedure  of  these  experiments  are:  1- 
evacuate  the  vessel  and  the  pipes  by  a  vacuum  pump  to  pressures  lower  than  80  millitorr; 
2-  assure  that  the  temperature  is  uniform  everywhere  and  fill  the  vessel  with  the  liquid  fuel. 
Then  close  the  vessel  valve.  3-  Evacuate  the  pipes  from  liquid  fuel;  4-  fill  the  vessel  with 
air;  5-  wait  at  least  for  three  minutes  to  let  the  mixture  becomes  homogeneous  and  uniform; 
6-  fire  the  vessel  with  ignition  system  and  record  the  pressure;  7-  vent  the  products  to  the 
atmosphere. 

Results:  The  autoignition  process  in  an  unburned  gas  mixture  is  a  very  fast  process 
with  intense  fluctuations.  Abnormal  pressure  fluctuations  can  be  a  trace  of  autoignition  in 
the  unburned  gas  zone.  One  of  the  important  issues  about  autoignition  is  its  initiation 
process.  In  an  ideally  homogeneous  mixture,  it  can  be  assumed  that  autoignition  happens 
everywhere  instantaneously  with  a  specific  ignition  delay  time  corresponding  to  the 
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temperature,  pressure,  and  equivalence  ratio  of  the  mixture.  In  these  conditions,  it  is 
assumed  that  the  mixture  is  perfectly  uniform  and  there  is  no  temperature,  pressure  or 
equivalence  gradient  in  the  mixture.  In  the  practical  combustors  like  internal  combustion 
engines,  it  is  hard  to  have  homogeneous  mixtures  since  there  is  a  certain  amount  of 
temperature,  pressure,  and  equivalence  ratio  gradients. 

In  this  research,  there  are  two  different  criterions  for  detecting  autoignition.  The 
first  is  done  using  the  pressure  signals  and  their  rates.  If  an  abrupt  jump  is  observed  in 
pressure  data  along  with  audible  noise,  it  can  be  a  trace  of  autoignition.  The  other  method 
is  by  analyzing  the  ionization  probes  signal.  If  an  ionization  probe  sends  out  highly 
oscillating  signals  (compared  with  nonnal  combustion)  it  will  be  a  sign  of  autoignition. 

Figure  21a  shows  normal  combustion  and  autoignition  of  JP-8/air  mixture  at  two 
different  initial  conditions.  It  is  seen  that  at  nonnal  combustion  the  pressure  increases 
smoothly  until  it  touches  the  vessel  wall,  reaches  to  equilibrium  condition  and  then  it 
decreases  because  of  energy  loss  to  the  wall.  The  peak  point  of  pressure  is  very  close  to  the 
amounts  calculated  by  the  STANJAN  equilibrium  code.  The  ionization  probes  signals  are 
nonnal  and  the  discontinuity  occurs  at  the  same  time,  so  the  flame  is  symmetrical  without 
any  buoyancy.  The  other  curve,  presents  the  autoignition  in  which  pressure  has  intense 
oscillations  along  with  a  very  high  peak  pressure  which  is  much  higher  than  corresponding 
equilibrium  pressure.  In  addition,  it  is  observed  that  the  ionization  probes  drive  out  very 
fluctuating  signals,  which  indicate  that  the  autoignition  has  occurred  in  the  end  gas  zone. 
Intense  fluctuations  in  the  ionization  probes  signals  are  due  to  autoignition  in  unbumed  gas 
zone  and  excitation  of  probes  before  flame  arrival.  Figure  21.b  shows  that  the  ionization 
signals  occur  almost  at  the  same  time.  Therefore  it  indicates  simultaneous  autoignition  all 
over  the  unburned  gas  zone.  After  repeating  the  experiments  and  ensuring  that  for 
reproducibility,  the  corresponding  pressure  and  temperature  should  be  measured  to  get  the 
critical  values.  Therefore  the  point  of  the  autoignition  should  be  found  to  measure  the 
corresponding  temperature  and  pressure.  This  can  be  done  with  two  different  approaches. 
The  first  method  is  by  analyzing  the  rate  of  pressure  variations  ( dP/dt )  values.  If  there  is 
any  abnormal  variation  in  dP/dt  values  compared  with  nonnal  combustion  dP/dt’ s,  it  can 
be  autoignition  trace.  Figure  22  is  the  corresponding  dP/dt  curve  of  Figure  21b.  It  is 
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observed  that  dP/dt  abrupt  increase  starts  almost  at  t  =  44.84  ms.  At  this  time  pressure 
increases  along  with  oscillations  that  can  be  caused  by  acoustic  pressure  waves. 

Figure  22  also  shows  the  ionization  probes.  It  indicates  that  the  signals  do  not  occur 
exactly  at  the  same  time  that  dP/dt  starts  to  increase  abnormally.  The  difference  is  about 
0.4  ms.  this  event  is  reproducible  in  most  cases.  It  can  be  explained  due  to  the  fact  that  in 
low  temperatures  the  ignition  process  is  in  two  stages  [14-16].  In  the  first  stage  of  ignition 
the  pressure  increases  slightly.  Then  in  the  second  part  of  the  ignition  where  the  main  part 
of  energy  is  released,  pressure  increases  at  a  higher  rate.  The  duration  of  these  stages 
depends  on  the  conditions  of  the  mixtures,  but  in  most  of  the  conditions  its  order  is  in 
milliseconds  [14-16]. 

As  it  is  seen  in  Figure  23,  the  difference  between  the  first  peak  of  pressure 
and  the  other  peaks  are  about  0.2-1  ms.  Therefore  ionization  probes  send  out  signals  when 
there  are  sufficient  radicals  and  products  in  the  mixture  to  ionize  the  wires  properly.  Since 
the  beginning  of  the  ignition  is  from  the  start  of  the  first  stage,  the  autoignition 
temperatures  and  pressures  are  measured  in  that  time. 

Figure  24  shows  the  temperatures  and  pressures  in  which  autoignition  took  place 
for  three  different  equivalence  ratios.  The  explosion  temperature  in  a  flame  propagating 
combustor  is  the  temperature  where  the  ignition  delay  time  of  the  unbumed  gas  mixture  is 
smaller  than  the  time  scale  of  flame  propagation. 

It  can  be  seen  that  there  is  a  specific  range  of  temperature  for  the  explosion  of  JP- 
8/air  mixtures.  It  varies  from  660K-710  K  in  a  wide  range  of  pressures  and  equivalence 
ratios.  This  figure  shows  the  sensitivity  of  the  mixture  to  temperature  is  much  more  than  to 
pressure.  The  range  of  autoignition  pressure  is  about  20  atm,  but  the  range  of  temperature  is 
just  about  40  K.  So  it  can  be  concluded  that  in  an  engine  where  JP-8  is  used  as  the  fuel  the 
temperature  does  not  exceed  the  temperature  range  presented  in  Figure  24. 
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Figure  la.  (top)  Schematic  of  15.24  cm  (6  inch)  ID  spherical  combustion  chamber 
and  lb.  (bottom)  13.33  cm  (5.25  inch)  ID  cylindrical  combustion  chamber,  with 
aspect  ratio  of  1 . 
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Figure  2.  Three  different  regions  of  gases  in  combustion  vessel 


P/Pi=4.97  '  P/Pi=3.08  P/Pi=2.32 


Figure  3:  JP8/air  mixture  0=0.8,  Ti=463  K,  a):  P;=l  atm,  b):  P;=2.5  atm,  c):  P;=5  atm 
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a  be 


P/Pi=2.34  P/Pi=2.26  P/Pi=2.97 

Figure  4:  JP8/air  mixture  0=0.9,  Ti=463  K,  a):  P;=l  atm,  b):  Pi=2.5  atm,  c):  P;=5  atm 
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P/Pi=2.62  P/Pi=2.4  P/Pi=2.14 

Figure  5:  JP8-air  mixture  0=1,  Ti=463  K,  a):  Pi=l  atm,  b):  P,=2.5  atm,  c):  P;=5  atm 
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P/Pi=2.6 


P/Pi=2.72 


P/Pi=l  .35 


Figure  6:  JP8-air  mixture  0=1.2,  Ti=463  K,  a):  P,=  l  atm,  b):  P;=2.5  atm,  c):  P;=5  atm 
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Pressure  (psi) 


Figure  7:  The  pressure  data  versus  time  for 
JP8-air  mixture,  P;=l  atm,  T ,=463  K 
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Volts 


TempemtuieiK) 


Figure  8:  JP8-air  mixture  laminar  burning  speed  at  different 
equivalence  ratios,  T ,=463  K,  Pj=l  atm 


Figure  9:  JP8-air  mixture  laminar  burning  speed  at  different 
equivalence  ratios,  T,=463  K,  P,=2  atm 


Figure  10:  JP8-air  mixture  laminar  burning  speed  at  different 
equivalence  ratios,  T ,=463  K,  P,=5  atm 
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a 


b 


P/Pi=2.47  P/Pi=2.52 

Figure  11:  JP8-09-diluent  mixture  0=1,  T ,=463  K,  Pi=2.5  atm,  a):  with  N2,  b):  with 


Argon 
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P/Pi=2 


P/Pi=l  .86 


P/Pi=l  .95 


P/Pi=2.52  1  P/Pi=2.68  1  P/Pi=3.46 

Figure  12:  JP8-02-He  mixtures,  T;=463  K,  P;=4  atm  a):  0=0.8,  b):  0=0.9,  c):  0=1 
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Figure  13:  JP8-C>2-He  mixture  laminar  burning  speed,  T;=463  K,  P;=l  atm 
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Figure  14:  JP8-C>2-He  mixture  laminar  burning  speed,  T;=463  K,  P;=2  atm 
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Figure  15:  JP8-C>2-He  mixture  laminar  burning  speed,  T;=463  K,  P;=4  atm 
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Figure  16:  JP8-C>2-He  mixture  laminar  burning  speed,  T,=463  K,  P;=5  atm 
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Figure  17:  JP8-C>2-He  mixture  laminar  burning  speed,  T;=463  K,  P;=7  atm 
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a 


P/Pi=2.47  P/P;=2.52  P/Pi=3.04 

Figure  18:  JP8-09-diluent  mixture  <3>  =  1,  Ti=463  K,  P,=2.5  atm,  a)  with  N2,  b)  with 

Argon,  c)  with  Helium 
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Figure  19:  Diesel/air  mixture  (p=  1,  T ,=463  K,  P;=l  atm 
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Figure  20:  Burning  speed  of  diesel/air  mixtures,  Pi=l  atm,  Ti=450  K 


0  10  20  30  40  50  60 


Time  (ms) 

(a) 


Time  (ms) 


(b) 

Figure  21:  Pressure  and  ionization  probes  data  of  JP-8/air,  O  =  1,  T;  =  483  K,  a:  P;=  1  atm, 

b:  Pj  =  10  atm 


25000 


-20000  Time  (ms) 

Figure  22:  Pressure  rate  JP-8/air,  O  =  1,  P;=  lOatm,  T;=  483  K 


Time  (ms) 

Figure  23:  Focused  pressure  of  JP-8/air,  O  =  1,  P;=  lOatm, 

T,=  483  K 
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750 


Pressure  (atm) 

Figure  24:  Explosion  temperature  versus  pressure  of  JP-8/air  mixtures 
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Table  1 :  Laminar  and  non-laminar  conditions  of  JP8-air  flames 

(p  Pi=l  atm  Pi=2.5  atm  Pi=5  atm 

Pcr=16.2  atm 

0.8  Non-Cellular  Non-Cellular 

Tcr=643  K 

Pcr=  10.2  atm  Pcr=  11.2  atm 

0.9  Non-Cellular 

Tcr=672  K  TCr=565  K 

Pcr=7.48  atm  Pcr=7.89  atm 

1  Non-Cellular _ 

Tcr=621  K  Tcr=504K 

Pcr=2.78  Pcr=3.87  atm 

1.2 


Tcr=594  K 


Tcr=541  K 


From  the 
beginning 


Table  3:  Cell  formation  analysis  of  JP8-C>2-He  mixtures 


<p 

Pi=2.5  atm 

Pi=4  atm 

0.8 

Pcr=  16.05  atm 

Pcr=  19.31  atm 

Tcr=832  K 

Tcr=760  K 

0.9 

P Cr=  15.17  atm 

Pcr=  15.71  atm 

Tcr=817  K 

Tcr=716K 

1 

Pcr=13.8  atm 

Pcr=13.3  atm 

Tcr=797  K 

Tcr=62 1  K 

43 


Table  5:  Cell  formation  and  Instability  Analyses  for  JP8-air  mixtures 


Table  6:  Cell  formation  and  Instability  Analyses  for  JP8-C>2-He  mixtures 


<p 

Pi  =  2.5  atm 

Pi  =  4  atm 

T  =  832  K 

T  =  760  K 

0.8 

P  =  16.05  atm 

P  =  19.34  atm 

a  =  0.0000395  m2/s 

a  =  0.0000278  m2/s 

5  =  0.0124mm 

8  =  0.01mm 

T  =  81 7  K 

T  =  716  K 

0.9 

P  =  223.855  psi 

P  =  231  psi 

a  =  0.0000405  m2/s 

a  =  0.0000312  m2/s 

8  =  0.01  mm 

8  =  0.01  mm 

T=797  K 

T=665  K 

1 

P=203  psi 

P=1 81  psi 

a  =0.000043  m2/s 

a  =.000035  m2/s 

8  =  0.01  mm 

8=0.0125  mm 
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